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ABSTRACT 

O 
O 

^ Sagittarius (Sgr) A* is a compact radio source at the Galactic center, pow- 

p I ered by accretion of fully ionized plasmas into a supermassive black hole of 

~ (3 — 4) X IO^Mq. However, the radio emission cannot be produced through 



X 
6 



the thermal synchrotron process by a gravitationally bounded flow (Liu & Melia 



2001). General relativistic magneto-hydrodynamical (GRMHD) simulations of 
black hole accretion show that there are strong unbounded outflows along the 
accretion. With the flow structure around the black hole given by GRMHD simu- 
lations, we investigate whether thermal synchrotron emission from these outflows 
may account for the observed radio emission and discuss the implications of this 
study on these GRMHD simulations and possible production of non-thermal par- 



c§ tides by this source. We find that simulations producing relatively high values 

of plasma j3 cannot produce the radio flux level without exceeding the X-ray up- 
^ per limit set by Chandra observations through the bremsstrahlung process. The 

00 predicted radio spectrum is also harder than the observed spectrum both for the 

^ one temperature thermal model and a simple nonthermal model with a single 

O power-law electron distribution. Since higher frequency emission is produced at 

smaller radii, the electron temperature needs to be lower than the gas temper- 
ature near the black hole to reproduce the observed radio spectrum. A more 
complete modeling of the radiation processes, including the general relativistic 
. ^ effects and transfer of polarized radiation, will give more quantitative constraints 

^ on physical processes in Sgr A* with the current multi-wavelength, multi-epoch, 

and polarimetric observations of this source. 
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Introduction 



Sagittarius (Sgr) A*, the compact radio source at the Galactic center, is powered by 
accretion of a supermassive black hole of ~ (3 — 4) x IO^Mq (Schodel et al. 



al. 2004) in the prevailing stellar winds in the Galactic center region (Melia 1992). The 



2002 



Ghez et 



bolometric luminosity of Sgr A* is more than 9 orders of magnitude lower than the corre- 



sponding Eddington luminosity, suggesting a radiatively inefficient accretion flow (Melia et 
al. [2001 Yuan et al. [2003 ). The radiative cooling processes therefore may be ignored while 
studying the dynamics of the accretion flow near the black hole, which simplifles the dy- 
namical equations and the corresponding numerical algorithms signiflcantly. Several general 
relativistic magneto-hydrodynamical (GRMHD) codes have been developed over the past 
few years to study the structure of non-radiative accretion flows near black holes quantita- 
tively (De Villiers et al. 2003; Gammie et al. 2003 Anninos et al. 2005). Given the 



extensive observations available for this supermassive black hole, it provides perhaps the 
best opportunity to study the physical processes in the strong gravity near black holes with 



GRMHD simulations (Falcke et al. 2000 Huang et al. 2008). 



The observed linear polarization and variability of the millimeter to near-infrared (NIR) 
emissions suggest that they are produced by the synchrotron process near the black hole 



(Aitken et al. 2000 Meha et al. 2000 Genzel et al. 2003; Eckart et al. 2004 2006). The 



longer wavelength emissions are circularly polarized (Bower et al. 2001) with weak linear 



polarization observed during frequent outbursts or flares (Zhao et al. 2004 Yusef-Zadeh 



et al. 2007), which in combination with the continuity of the centimeter to sub- millimeter 



spectrum suggests a synchrotron origin for the longer wavelength emissions as well (Falcke 



et al. 1998 An et al. 2005). The analogy of Sgr A* with radio loud AGNs also favors a 
synchrotron scenario for the radio emission. However, the longer wavelength flux densities 
are less variable than the millimeter and sub-millimeter flux densities, indicating emissions 



from relatively larger radii (Falcke & Markoff 2000; Liu & Melia 2001 Herrnstein et al. 



2004). Indeed, the millimeter and higher frequency emissions have been explained reasonably 
well with an accretion torus within ~ 10 rs of the black hole, where rs is the Schwarzschild 



Goldston et al. 


2005 


Liu et al. 


2007 


Noble et al. 


2007 



high spatial resolution VLBI observations do show that the intrinsic size of the millimeter 
emission region is within 20 rs in radius (Bower et al. 2004; Shen et al. 2005). While 



3 



the general relativistic (GR) effects are important for these shorter wavelength emissions 
originating near the black hole (Falcke et al. 2000; Bromley et al. 2001), the longer 



wavelength emissions are produced beyond ~ 10 so that the GR effects are insignificant 
and one may readily use the GRMHD simulations to model the observed emission without 
using the sophisticated GR ray-tracing radiation transfer codes for the polarized synchrotron 



emission (Broderick & Loeb 2005 Huang et al. 2008). 



High resolution X-ray observations with the Chandra space telescope shows that the 
quiescent X-ray excess from the direction of Sgr A* is extended and has ion emission lines 
in the spectrum, suggesting that the emission is mostly produced at large radii near the 



capture radius of the accretion flow by a plasma of a few keV in temperature (Baganoff et 



al. 2001 2003 Xu et al. 2006). X-ray flares are routinely observed from the direction 



of Sgr A* (Belanger et al. 2006). Their spectra can be fitted with a featureless power law 



(Porquet et al. 2003), and they appear to be always accompanied by NIR flares, which are 



produced near the black hole (Eckartetal. 2004 2006 Yusef-Zadeh et al. 2006 2008). The 



correlation between the NIR and X-ray flares indicates a synchrotron self-Comptonization 
(SSC) origin for the X-ray flares ( [Markoff et al. ||2001| pu fc Meha |[2001| ), though the X-ray 
flares may also be produced via the bremsstrahlung process enhanced by instabilities related 



to the accretion process (Liu & Melia 2002 Tagger & Melia 2006). Recent observations of 



flares seem to favor a synchrotron process (Liu & Melia 2001; Dodds-Eden et al. 2009). 



The quiescent state X-ray emission from the inner accretion flow therefore must be lower 
than the observed flux level. X-ray emission can be produced through both the SSC and 
bremsstrahlung processes. The observed high radio flux level and low upper limit of X-ray 
flux suggest that the radio emission cannot be produced via thermal synchrotron emission in 
a bounded flow ( Liu fc Melia ||200T ). It therefore has to be produced by unbounded outflows 
and/or by nonthermal populations of relativistic electrons. Given the non-radiative nature 
of the GRMHD simulations, strong outflows and winds appear to be inevitable, especially 



for highly spinning black holes (Noble et al. 2006 Hawley et al. 2007). In this paper, we 



investigate whether emissions from these simulated accretion flows can count for observations 
of Sgr A*, especially the low frequency radio emissions, self-consistently. 

We carried out 2-dimensional (2D) GRMHD simulations with the Cosmos++ and 



HARM codes (Gammie et al. 2003 Anninos et al. 2005) assuming axis-symmetry of 



the accretion flow. The accretion disk therefore must be aligned with the equatorial planes 
of spinning black holes. The length scale is determined by the black hole mass. Since no 
radiative cooling is included, the density is scaled linearly with the mass accretion rate, and 
the gas temperature and plasma /3, deflned as the ratio of the gas pressure to the mag- 
netic field pressure, are independent of the black hole mass and accretion rate. GRMHD 



-4- 



simulations give the gas density, pressure, and magnetic field. The gas temperature can 
be derived from the equation of state. The mass of the supermassive black hole in Sgr A* 



is well-measured by observations of stellar orbits around the black hole (Schodel et al. 



2002 Ghez et al. 2004). One therefore needs to adjust the accretion rate, black hole spin. 



and electron distribution to fit the observed spectrum of Sgr A*. We consider both a one 
temperature and a two-temperature model of electrons in the thermal synchrotron emission 
scenario. For non-thermal synchrotron models, the electron distribution is assumed to follow 
a steep power-law with a low energy cutoff. In § 2 we discuss the methods that we used to 
model the emission spectra from the simulated accretion flows. A detailed discussion of our 
results and the implications are given in § 3. 



2. Modeling the Emission Spectrum 



We run the Cosmos++ code over a simulation domain of 218 in radius. Initially, 
there is an equilibrium torus with the pressure maximum at 100 rs- Under the influence of 
the strong gravity of the central black hole, the torus evolves and the magneto-rotational- 
instability sets in and induces the accretion process (Balbus & Hawley 1991). Specific 



attributes of the accretion disc were examined and the magnetic field, gas temperature, 
and particle density profiles are shown in Figure [TJ Because the simulation was run up 
to the dynamical time near 100 rs, the initial torus structure has not relaxed completely 
and, beside the cusped accretion torus in the inner region, which corresponds to the particle 
concentration within 100 rs, there is another high density region at ~ 100 rs, which maps 
to the initial torus. If the simulation were allowed to run longer, this region of high density 
would be pulled closer to the black hole and merge with the inner accretion torus. The 
simulation also produces extremely high temperature regions along the polar directions, the 
polar regions within a half open angle of 6° and the inner 10 rs region are excluded in our 
modeling of the emission for this simulation. 

We first consider the simplest case, where electrons and ions reach thermal equilibrium 
throughout the simulation domain, so that the electron temperature is given by the gas 



temperature. The synchrotron emission coefficient at frequency u is given by (Pacholczyk 



1970] [Mahadevan et al. ][I996| |Liu et al. |[2006| ) 



where 
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Fig. 1. — Left: The magnetic field in units of Gauss in a logaritlimic color scale for the 
simulation with the Cosmos++. Middle: Same as the left panel but for the temperature 
in Kelvin. Right: Same as the left panel but for the particle number density in cm~^. 
The outermost high density region corresponds to the position of the initial torus. If the 
simulation were to run longer this region of high density at about 100 would merge with 
the inner high density region. 

I{xm) = 4.0505x-^/^(l + 0.40x-^/^ + O.SSlGx^^/^) x exp(-1.8899x^/^) , (3) 
1c = ^, (4) 

where = T = PgasTnp/2pk^ is the electron temperature and T, Pgas, P, B, n, e, rrie, c, 
and /cb indicate the gas temperature, pressure, mass density, magnetic field, particle density, 
the elementary charge unit, electron mass, speed of light, and the Boltzmann constant, 
respectively. The left and middle panels of Figure |2] show the profiles of £^ at 1.0 and 100 
GHz, respectively. The corresponding absorption coefficient hz^ and the opacity r along the 
light of sight of the accretion disk is given, respectively, by 

= „ y , T = [ K^dl, (5) 

where the integration of r is along the light of sight from the observer into the source region. 



Then the observed flux density from the accretion flow (|Melia et al. 2001|) 

L 



F.{iy) = ^ I I^ds, (6) 



Fig. 2. — Left: The synchrotron emission coefficient in the plasma co- moving frame at 
1.0 GHz in the c.g.s. units. Middle: Same as the left panel at 100 GHz. Right: The 
bremsstrahlung emission coefficient at lO^'^Hz, which corresponds to ~ 4.1 keV in energy. 



where 



J £y exp(— r) dl 



is the specific intensity and D is the distance to the Galactic Center and the integration of 
is over the projected area of the source. For a face-on disk with the observer located 
aX z = the radio spectra are shown in the left panel of Figure |3] for several values of 
the density normalization. These values are chosen so that the model spectra embrace the 
observed radio flux densities. The red line with a 10 times higher density proflle than that 
shown in Figure [T] gives the best fit to the observed spectrum. 

However, the high gas density implies strong X-rays via the bremsstrahlung process 
from the two torii evident in the density profile of Figure [Tj The bremsstrahlung emission 
coefficient is calculated by (Rybicki & Lightman 1979[ ): 



£^ = 8.5 X 10-^VT-^/^ exp{-hu/kBTe) , (7) 

where h is the Planck constant. The right panel of Figure [2] shows the bremsstrahlung 
emission coefficient at ~ 4.1 keV and the corresponding emission spectra are also indicated 
in the left panel of Figure |3} The bremsstrahlung emission exceeds the observed upper limit 
for all the three density normalizations. This simple model therefore cannot produce the 
observed radio flux level without exceeding the X-ray upper limit. 
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It is well-known that nonthermal synchrotron emission is more efficient than thermal 
synchrotron. We next study whether a simple nonthermal model with a single power-law 
electron distribution can account for the radio flux without violating the X-ray upper limit. 
The electron distribution is given by 

N{E) = NqE-p for E> Eo, (8) 

where Eq is the low energy cutoff. A^^o and Eq can be calculated with 

NoE-PdE = n, (9) 

POO 

/ NoE~P+^ dE = Snk^T^ . (10) 

JEn 



Eo 



The corresponding emission and absorption coefficients are given respectively by (Pacholczyk 
[l970| ) 

£. = c,ip)No[Bsma]^P^'y'l^—j , (11) 

/ \ -(P+4)/2 

K, = c,ip)No[Bsma]^P^'^/'l^—\ , (12) 

where a is the angle between the magnetic field and the line of sight and c, are defined in 



Pacholczyk (1970). In our calculations we chose p = 3 to avoid significant NIR emission in 
the quiescent-state. As the p increases, the non-thermal emission model became less efficient 
and closer to the thermal model. The right panel of Figure [3] show the spectra for several 
density normalizations, where the synchrotron spectrum is cut off artificially at lO^^Hz and 
the bremsstrahlung spectra is given by the thermal formula of equation ([T]). Compared 
with the thermal model, more radio emission is produced for a given density normalization. 
However, to produce the observed radio flux level, the bremsstrahlung emission still exceeds 
the observed upper limit. 

While the outer torus of the simulation is expected to merge with the inner torus if the 
simulation is run for a long enough time, the inner torus is already relaxed at the present 
epoch of the simulation. Therefore a longer run may remove the bremsstrahlung emission 
from the outer torus, we don't expect the overall bremsstrahlung X-ray flux to decrease by 
more than a factor of 2. It is therefore challenging to produce the observed radio flux level 
from Sgr A* without exceeding the X-ray upper limit with this simulation. The high X-ray 
flux as compared with the radio flux is a direct consequence of the high values of plasma /3 
given by this simulation. 
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For this simulation, f3 is always greater than 100 except in the polar regions, where the 
GRMHD simulation is not reliable and we have excluded these regions while obtaining the 
emission spectra. The bremsstrahlung emission is determined by the gas temperature and 
density. The synchrotron emission also depends on the magnetic field. For high values of 
P, the magnetic field is weak. It is therefore difficult to produce significant radio emission 
without making X-rays through the bremsstrahlung. Simulations that are capable of pro- 
ducing lower values of P are needed to overcome this challenge. We also noticed that the 
radio spectra of both the thermal and non-thermal model are too hard to fit the observed 
spectrum. Decreasing the value of plasma (3 alone may not explain the observations. 

Errors in the divergence of B are handled differently in the Cosmos++ than the HARM 



code (Gammie et al. 2003 Anninos et al. 2005). The latter uses a contrained transport 



schemes for evolving the induction equation whereas Cosmos++ does not. As a consequence, 
the HARM code appears to be able to produce lower values of plasma f3. We ran the HARM 
code in a domain of 100 rs in radius with the torus pressure maximum initially at ~ 6 r5. 
The simulated disk structure in a late time when the initial torus has completely relaxed is 
shown in Figure |4} The simulation produces more than 10 times stronger magnetic fields 
with comparable values for the density maximum. Figure |5] gives the corresponding emission 
coefficient maps. Although the synchrotron emission may be comparable to that in Figure |4| 
we expect much less bremsstrahlung X-rays due to the compactness of the torus in this latter 
simulation. The left panel of Figure |6] shows the spectra of the thermal model with several 
density normalizations. The bremsstrahlung X-ray emission is far below the observed upper 
limit. However, the radio spectra are still not fiat enough to fit the observed spectrum. 

Higher frequency emission is produced at relatively smaller radii, where the magnetic 
field is strong and the gas temperature is high. The fact that the model predicted spectrum 
is harder than the observed spectrum suggests that emission from small radii needs to be 
suppressed to produce a softer spectrum. The synchrotron cooling time of electrons at 
smaller radii is shorter. It is possible that the electron temperature is actually lower than 
the gas temperature at small radii. Modeling of the emission from a small accretion torus 



also favors a lower electron temperature than the gas temperature (Liu et al. 2007 Noble 



et al. 2007 Huang et al. 2008). We model the electron temperature with 



(14) 



where R is the radius, beyond which electrons reach thermal equilibrium with the ions and 
5 is a free parameter determined by the electron heating and cooling processes. 



There was no need to apply the two-temperature model to the first simulation be- 
cause it would decrease the amount of NIR and radio emissions without decreasing the 
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bremsstrahlung emission values, we would still have the same problem of exceeding the up- 
per limits imposed by Chandra on the X-ray emission. We applied this two-temperature 
model with R = 10 rs and S = 1.0 to the second simulation and calculated the thermal 
emission spectrum. The right panel of Figure [6] shows the results. The model with a 1000 
times higher density than that in Figure |4] fits the radio spectrum reasonably well. The 
underproduction of radio emission below 10 GHz is due to the relatively small simulation 
domain of 100 rs- Most of the longer wavelength emission is produced beyond 100 rs- How- 
ever, the model predicted X-ray flux level just matches the X-ray upper limit. So this model 
flts the observations marginally. Including the Doppler effect and considering the possibility 
of nonthermal emission may improve the fltting. 



3. Conclusions 

Sgr A* is a radio source powered by accretion into a black hole. Modeling the emission 
from this radio source at the Galactic center provides us with knowledge of the dynamics of 
the accretion flow near black holes and helps us further understand the plasma physics pro- 
cesses of the accretion torus. We were unsuccessful in our attempt to reproduce the observed 
spectrum using simulations with very high values of plasma (3- With the HARM code, lower 
values of plasma P are produced. Although we were able to apply a two temperature model 
to this simulation to correctly match the shape of the observed radio emission spectrum, 
we had the reoccurring problem of exceeding the upper-bound imposed by Chandra in the 
X-ray emission spectrum. We also note that the HARM simulation has a much smaller 
initially torus, which effectively reduces the bremsstrahlung X-ray flux. Although this itself 



may suggest that the angular momentum of the accretion flow must be low (Coker & Melia 



1999 Rockefeller et al. 2004), one needs to explore the dependence of the bremsstrahlung 



X-ray emission on the size of the initial torus to quantify this constraint. 

The above challenge of modeling observations of Sgr A* with GRMHD simulations may 
be alleviated by several means. First, as we have shown in the paper, nonthermal synchrotron 
emission is more efficient than thermal synchrotron and in the nonthermal scenario the radio 
spectrum may be fitted with a lower density normalization than the thermal one, which 
reduces the bremsstrahlung X-ray fiux. The Doppler effect, which we didn't consider in the 
current study, is expected to boost the radio emission as well. Introducing a black hole spin, 
one can produce even lower values of plasma (3 and higher Doppler boosting, both of which 
will enhance the radio emission. Lastly, the bremsstrahlung X-ray may be suppressed by 
reducing the mean density in the accretion disk through the tilt of the accretion disk with 
respect to the black hole spin ( Liu fc Melia |[2002b ) . Three dimensional simulations are then 



needed to quantify this effect (Rockefeller et al. 



2005 



Fragile et al. 



2007) 



In our modeling, we haven't studied the dependence of the model spectrum on the 
inclinational angle of the accretion disk, the GR effect, and the transfer of polarized emission. 
However, the radio spectrum is expected to have very weak dependence on these factors. 
Nevertheless, they should be included in a more complete study. 
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Fig. 3. — Left: Thermal synchrotron and bremsstrahlung emission spectra for the one 
temperature model. The black, red, and blue lines correspond to the particle density in 
Figure [1} 10, and 100 times this density, respectively. As expected, the bremsstrahlung flux 
level scales with the density squared. The radio to sub-millimeter data was obtained from 
observations over the past decade. The scattering of the millimeter and sub-millimeter data is 
caused by the variability of the source. The NIR flux densities are for flares observed during 
the past few years and therefore should be considered as upper limit for the "quiescent" 
state emission from the accretion flow. The butterfly is given by Chandra observations of 
the quiescent state X-ray emission, which appears to be dominated by thermal emission at 
large radii. Thus it also should be interpreted as an upper limit for X-ray emission from the 
small scale accretion flow. All the bremsstrahlung spectra exceed this upper limit for X-ray 
emission. Right: Same as the left panel but for the non-thermal model with a single power- 
law electron distribution with a low energy cutoff. The synchrotron spectra are showed up 
to 10^^ Hz. The density normalization is indicated in the figure. Compare with the thermal 
models in the left panel, the nonthermal model is more efficient in the synchrotron emission 
process. However, to match the observed radio flux level, the bremsstrahlung X-ray flux still 
exceeds the observed upper limit. 
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Fig. 5.— 



Same as Figure |2] but for the simulation of Figure |4j 
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Fig. 6. — Left: The same as the left panel of Figure [3] but for the simulation of Figure |4j 
Right: The same as the left panel except that electron temperature is lower than the gas 
temperature toward small radii as given by equation (14) with R = 10 rs and 6 = 1.0. 



